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We present a Monte Carlo study of the correlation length £ of randomly diluted antiferromag- 
netic Heisenberg ladders, composed of two spin-1/2 chains. For weak and intermediate inter- 
chain couplings, J±/J < 1, we find an enhancement of correlations that is strongest for a fraction 
z* ~ J±/ (8 J) of dilutants. We are able to access the experimentally relevant low-temperature 
regime, T/J ~ 1/500, and find that the recently inferred Neel temperature of Sr(Cui_ z Zn z )203 
corresponds to a curve of constant correlation length £ « 18 of the single diluted ladder with 
J±/J J5 1/2- The primary reason for the Neel ordering is argued to be a strong enhancement of 
two-dimensional correlations due to a Cu-Sr-Cu exchange coupling of » WmeV in the stacking 
direction of the ladders. 
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Low-dimensional quantum Heisenberg antiferromag- 
nets exhibit many unusual collective properties owing to 
the increased importance of quantum fluctuations when 
the dimensionality of the system is less than three, and to 
the continuous 0(3) spin rotational symmetry of Heisen- 
berg interactions. The one-dimensional (ID) Heisenberg 
chain and its two-dimensional (2D) analog, the square- 
lattice antiferromagnet, have been studied extensively for 
many decades. Arguably, the most prominent examples 
in this context are the lamellar copper oxides, which in 
their undoped states are square-lattice Heisenberg anti- 
ferromagnets with predominant nearest-neighbor (NN) 
interactions and a spin value of 5 = 1/2 |Q. The prop- 
erties of the NN Heisenberg chain are very different from 
those of the square lattice. Instead of true long-range 
order, the 5 = 1/2 chain exhibits a critical ground state 
characterized by algebraic decay of the spin correlations 
and gapless excitations Q. In 1983, Haldane || conjec- 
tured that the basic properties of the 5 = 1/2 chain are 
shared only by chains with half-odd-integer spin, and 
that integer-5 chains should behave fundamentally dif- 
ferently. The latter were predicted to possess a spin- 
liquid ground state characterized by a finite correlation 
length £ due to the presence of a gap A ~ l/£ to 
the lowest excitations. These predictions have been con- 
firmed both numerically and experimentally Q. 

In recent years, the question of how quantum fluctua- 
tions manifest themselves in systems whose extent is fi- 
nite in one dimension and infinite in another has received 
much attention. In particular, arrays of coupled chains 
of width n, so-called spin ladders, have been studied ex- 
tensively H . Good physical realizations of these systems, 
with 5 = 1/2, were found to exist in materials such as 
SrCu 2 3 (n = 2) and Sr 2 Gu 3 05 (n = 3) @. In close 
analogy to the behavior of the spin-5 chains, the prop- 
erties of NN 5=1/2 ladders are fundamentally different 
for even and odd numbers of coupled chains || . Quite 
generally, it is believed that quantum fluctuations always 
generate a spin gap when the product nS is an integer, 



and that the destructive interference of fluctuations leads 
to a gapless spectrum for half-odd-integer values of nS 
. Theory for the temperature dependence of the spin 
correlations of 5 = 1/2 ladders ||] is in good agreement 
with Monte Carlo results J||| ■ 

A surprising recent experimental discovery has 
been that random dilution in the ladder material 
Sr(Cui_ z Zn z ) 2 3 (5 = 1/2, n = 2) leads to antiferro- 
magnetic long-range order at low temperatures, even for 
very small Zn concentrations (To). A schematic of the 
structure of this material is shown in Fig. 1(a). The 
observed Neel order is counter-intuitive, since based on 
classical intuition one might have expected a diminution 
of correlations due to dilution, and hence a reduced ten- 
dency to order. A large number of numerical and theo- 
retical studies have addressed this issue [jll| , yet a quanti- 
tative understanding of the experimentally relevant low- 
temperature spin correlations of a diluted two-chain lad- 
der is still lacking. Furthermore, the explicit pathway 
by which dilution promotes three-dimensional (3D) Neel 
order remains elusive. 

In this paper, we present a numerical investigation 
of the effects of random dilution z on the spin corre- 
lations of a n — 2, S = 1/2, NN antiferromagnetic 
Heisenberg ladder. We are able to reach the experimen- 
tally relevant low-temperature regime, T/J ks 1/500, 
and find that a small amount of dilution indeed en- 
hances the antiferromagnetic correlation length, albeit 
rather weakly. For weak and intermediate inter-chain 
couplings, Jj_/J < 1, we find that the enhancement is 
strongest for an average spacing between dilutants that 
approximately equals £o = £,( z = 0, J_i_/ J, T/J = 0), the 
zero-temperature correlation length of the un-diluted 
ladder. In the limit of small J±/J, we have shown pre- 
viously that £o ~ 3.9J/Jj_ ||- Since the Neel temper- 
ature T N (z) of Sr(Cui_ z Zn z ) 2 3 peaks at z w 4% [@, 
this in turn requires that Jj_/J < 1/2. This bound is 
consistent with recent experimental and theoretical esti- 
mates for the ratio J±/J. Moreover, we find that Tjv(z) 
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of Sr(Cui_ z Zn z )203 corresponds to a curve of constant 
correlation length £ « 18 of a single ladder. We argue 
that the primary reason for the experimentally observed 
divergence of spin correlations at Tm{z) is a crossover 
from ID to 2D behavior due to the inter-ladder Cu-Sr- 
Cu exchange in the stacking direction of the ladders; 
an important observation is the fact that the Cu-Sr- 
Cu geometry and distances are very similar to those in 
the Cu-Y-Cu bi-layer in the high-T c precursor material 
YBa2Cu306+5- From neutron scattering measurments 
(for S = 0.15 — 0.20) it is known that the bi-layer cou- 
pling is w lOmeV E|. The final 3D antiferromagnetic 
ordering then occurs because of a very weak pseudo- 
dipolar coupling between the frustrated ladder segments 
jy§. We thus predict a broad ID to 2D crossover of 
the spin correlations followed by sharp 3D Neel ordering 
transition. This could be observed directly in neutron 
scattering experiments once suitably large single crystals 
become available. 

As in our previous studies the ladders are inves- 
tigated with a very efficient loop cluster algorithm |fl4) . 
The Hamiltonian operator for a Heisenberg ladder is 



H 



(ij) 



(ij)± 
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where Si = \(Ji is the quantum spin operator at each site 
% of non-zero spin, while (ij) and (ij)x denote nearest 
neighbors along the chains and across the rungs of the 
ladder, respectively (see Fig. 1(b)). The couplings con- 
sidered are antiferromagnetic, that is, J and Jj_ > 0. Pe- 
riodic boundary conditions are employed only along the 
chains, and we use units in which the lattice constant is 
unity, as well as % = kg — g^B = 1. We typically average 
over 40 random configurations of each 10 4 updates and 
measurements, and compute the staggered instantaneous 
correlation function 
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where sign(i, j) = 1 if the spins at i and j are on the same 
sublattice of the antiferromagnet, and sign(i,j) = —1 
otherwise. The correlation length of the diluted lad- 
ders, £ = £(z, J±/ J,T/J), is extracted at large dis- 
tances r, r/£ > 3, from fits to the symmetrized form 
C(r) ~ cosh[(L/2 — the length L of the ladders is 

kept at least an order of magnitude larger than £. 

In the present study, we focus on ladders with 1/4 < 
J±/J < 1. The dependence of the correlation length 
on the degree of dilution is shown in Fig. 2. Here, the 
temperature is chosen such that T/J± = 3/100 is con- 
stant. As z is increased from zero, £ is enhanced, reaches 
a maximum, and eventually decreases at large z. In a 
previous study of un-diluted ladders || we found that 
at weak inter-chain coupling the product A£ exhibits a 
rather simple scaling behavior in the variable T/A. The 



weak-coupling regime is characterized by the relation- 
ship A w 0.41Jj_, and it sets in below J ± /J rj 1/2 
This motivated us to re-plot the data of Fig. 2 as £ J± / J 
versus zJ/Jj_, as shown in the inset. Indeed, the data for 
J± = 1/4 and 1/3 collapse to a single curve, while those 
for J± > 1/2 deviate slightly from this curve. 

The enhancement of correlations appears to be largest 
for an approximate average spacing between dilutants of 
£o = £(0, J±/J, 0); specifically we find empirically that 
the maximum occurs for z* — | £ a +£* i where £* = 0(1) 
is a constant. The weak-coupling regime is characterized 
by A w 0.41Ji and £ = cm/ A, where cm — 7i72 is the 
spin-wave velocity of the single spin-1/2 chain M. Deep 
in the weak-coupling regime, £o ~ 3.9J/ J± is large, so 
that z* w 0.13Jj_/J, consistent with Fig. 2. 

For Sr(Cui_ z Zn z )203, the Neel temperature, and thus 
the enhancement of correlations, is the largest for z as 
4%. Based on our observations for the single ladder, this 
implies that 1/4 < J±/J < 1/2. Such a strong exchange 
anisotropy might at first appear to be somewhat surpris- 
ing, given that the Cu-O-Cu distances along the chains 
and across the rungs are nearly identical. However, both 
recent experiments [|l6|l7 and theoretical calculations 
Q imply that J±/J is of order 0.4 to 0.5, consistent 
with our own deductions. 

We now discuss quantitatively the Neel ordering of 
Sr(Cui_ z Zn z ) 2 03, based on the assumption that J±/J w 
1/2. To effect this, we simulate the diluted ladder with 
Jj_/J = 1/2 down to very low temperatures, T/J w 
1/500. In our previous study of the the un-diluted ladder 
(z = 0) with J ± /J = 1/2 we had found that £ ~ 7.5 [|. 
The temperature dependences of the correlation length at 
several non-zero values of z are shown in Fig. 3 together 
with the previous result for z — 0. Corresponding stag- 
gered susceptibility data are given in the inset. As noted 
before, the enhancement of correlations is strongest for 
z* a 4%. 

The maximum ordering temperature of the experimen- 
tal system is T N a 8K at z* « 4% p| , as shown in 
Fig. 4. With J « 1900if @, this corresponds to 
Tn/J « 0.004. From Fig. 3, it is apparent that the corre- 
lations of the isolated ladder for z « 4% and T/J as 0.004 
are enhanced only to £ » 18. In fact, we find that Tm{z) 
of Sr(Cu z Zni_ z )203 corresponds to a curve of constant 
correlation length, £ = 18(1), as shown in Fig. 4. 

What then is the mechanism that leads to the dra- 
matic enhancement of correlations and the concomitant 
Neel order observed experimentally? The fact that the 
experimentally observed Tn(z) corresponds to a curve of 
constant correlation length of a single diluted ladder im- 
plies that this mechanism is independent of the degree of 
dilution. It is important to realize that the inter-ladder 
Heisenberg coupling within the Cu-0 sheets remains es- 
sentially frustrated at non-zero dilution. A key obser- 
vation in understanding the onset of long-range order 
in Zn-doped SrCu 2 3 is that the Cu 2+ - Sr 2 + - Cu 2+ 
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geometry (see Fig. 1(a)) and distances in the stack- 
ing direction of the ladders are essentially identical to 
those in the Cu 2+ — Y 3+ — Cu 2+ bi-layer arrangement 
in YBa2Cu306+a; recent neutron scattering experiments, 
for S = 0.15 - 0.2, give J c w lOmeV @. There- 
fore, the primary perturbation to the quasi-lD ladders 
in Sr(Cui_ z Zn z )203 is the coupling between ladders in 
adjacent planes. The reduced value for this inter-ladder 
coupling is ct2D = Jc/J ~ 6%. Accordingly, we antici- 
pate a progressive crossover from ID to 2D spin correla- 
tions, followed by a sharp transition to 3D Neel order due 
to a intra-planar pseudo-dipolar coupling between ad- 
jacent ladders Jl^|. This evolution should be observable 
experimentally once suitably large single crystals become 
available. 

The inset of Fig. 3 shows some results for the stag- 
gered susceptibility per spin, Xs ~ T^ 1 ((J2i{—^Y S*) 2 ) ■ 
We find that J\s ~ 160 along the curve of constant cor- 
relation length £ k, 18, Fig. 4. Mean-field theory for the 
enhancement of the susceptibility predicts 



which diverges for J\ s = 1/ck2D ~ 18. From the inset 
of Fig. 3 it is clear, that such a small value of J\s is 
already reached at rather high temperatures. Perhaps 
not surprisingly, mean-field theory predicts 2D ordering, 
and, implicitly, 3D ordering, at a temperature which is 
about an order of magnitude higher than that observed 
experimentally for all values of z. Empirically, the condi- 
tion for the Neel ordering appears to be £ = l/ain ~ 18, 
rather than J\ s = I/o^d; this remains to be understood 
theoretically. 

A previous numerical study of the susceptibilities of 
the diluted ladder |l!| extended to temperatures as low as 
T/J = 1/200 and used Eq. (3) together with J ± /J=l, 
J = lOOOif, and a^D ~ 2%, to explain qualitatively the 
phase diagam of Sr(Cui_ z Zn z )203. Our results are con- 
sistent with this work, but give a physically more com- 
plete and precise description. 

Recently, Fujiwara et al. [ ^of presented an interesting 
NMR study of 65 Cu linewidths in Sr(Cui_ z Zn z ) 2 3 at 
small Zn concentrations, z < 0.5%, and and for tempera- 
tures T > 20K. Our results imply that these data cannot 
be understood with the model suggested by the authors 
of Ref. ^0|, since the latter model requires ladder cor- 
relation lengths one to two orders of magnitude larger 
than those which we compute. Further work is required 
to understand these most interesting NMR results. 

In summary, we have numerically determined the cor- 
relation length of randomly diluted spin-1/2 two-chain 
Heisenberg ladders down to the experimentally relevant 
low-temperature regime, T/J « 1/500. At weak and in- 
termediate inter-chain couplings, Jj_/J < 1, we find an 
enhancement of correlations that is strongest for a frac- 
tion z* w J±/(8J) of dilutants. The recently inferred 



Neel ordering temperature Tn(z) of Sr(Cui_ z Zn z ) 2 03 is 
found to correspond well to a curve of constant corre- 
lation length £ » 18 of the single diluted ladder with 
J±/J 1/2. We argue that the primary reason for the 
experimentally observed 3D Neel ordering is an enhance- 
ment of 2D correlations due to an antiferromagnetic cou- 
pling OL2D — Jc/J ~ 6% in the stacking direction of the 
ladders. Empirically, the condition for Neel ordering thus 
appears to be £ = \/a.iD- The ID to 2D crossover is then 
followed by a sharp transition to 3D Neel order due to a 
tiny intra-planar pseudo-dipolar coupling between adja- 
cent ladders. Once sizable single crystals become avail- 
able, this predicted evolution of the correlations should 
be observable in neutron scattering experiments. 

We would like to thank A. Aharony, H. Fukuyama, T. 
Imai, D.C. Johnston, Y.-J. Kim, P.A. Lee, T.M. Rice, 
and U.-J. Wiese for stimulating discussions and corre- 
spondence. This work was supported by the NSF under 
Grant No. DMR 97-04532, the Joint Services Electronics 
Program, and the International Joint Research Program 
of NEDO (New Energy Industrial Technology Develop- 
ment Organization, Japan). 



[10] M. Azuma et al, Phys. Rev. B 55, R8658 (1997). 

[11] See for example: N. Nagaosa, A. Furusaki, M. Sigrist, 
and H. Fukuyama, J. Phys. Soc. Jpn. 65, 3724 (1996); Y. 
Motome, N. Katoh, N. Furukawa, and M. Imada, J. Phys. 
Soc. Jpn. 65, 1949 (1996); Y. lino and M. Imada, J. Phys. 
Soc. Jpn. 65, 3728 (1996); H.-J. Mikeska, U. Neugebauer, 
and U. Schollwock, Phys. Rev. B 55, 2955 (1997); G.B. 
Martins, M. Laukamp, J. Riera, and E. Dagotto, Phys. 
Rev. Lett. 78, 3563 (1997). 

[12] D. Reznik et al, Phys. Rev. B 53, R14741 (1996); S.M. 
Hayden et al, Phys. Rev. B 54, R6905 (1996). 

[13] The pseudo-dipolar coupling in the ladder materials 
is closely analogous to that in the frustrated square- 
lattice antiferromagnet S^CusdClj; see F.C. Chou et 
al, Phys. Rev. Lett. 78, 535 (1997). 



[1] M. Greven et al, Z. Phys. B 96, 465 (1995). 

[2] H.A. Bethe, Z. Phys. 71, 205 (1931); A. Luther and I. 

Peschel, Phys. Rev. B 9, 2911 (1974). 
[3] F.D.M. Haldane, Phys. Lett. A 93, 464 (1983). 
[4] See, e.g., I. Affleck, Rev. Math. Phys. 6, 887 (1994). 
[5] For a review, see E. Dagotto and T.M. Rice, Science 271, 

618 (1996). 

[6] Z. Hiroi et al, J. Solid State Chem. 95, 230 (1991). 

[7] D.V. Khveshchenko, Phys. Rev. B 50, 380 (1994); SR. 

White, R.M. Noack, and D.J. Scalapino, Phys. Rev. Lett. 

73, 886 (1994); G. Sierra, J. Phys. A 29, 3299 (1996). 
[8] O.F. Syljuasen, S. Chakravarty, and M. Greven, Phys. 

Rev. Lett. 78, 4115 (1997). 
[9] M. Greven, R.J. Birgeneau, and U.-J. Wiese, Phys. Rev. 

Lett. 77, 1865 (1996). 



3 



[14] H.G. Evertz, G. Lana, and M. Marcu, Phys. Rev. Lett. 
70, 875 (1993); U.-J. Wiese and H.-P. Ying, Z. Phys. B 
93, 147 Q-994). 

[15] In Ref. M, we obtained A = 0.41 (1)J±; based on addi- 
tional data for J± — J/3 and J/4 in the present study, 
we find A = 0.407(5) J±. 

[16] D.C. Johnston, Phys. Rev. B 54, 13009 (1996), and pri- 
vate communication. 

[17] R.S. Eccleston et al. (preprint). 

[18] M. Arai and H. Tsunetsugu, Phys. Rev. B 56, R4305 
(1997); T.F.A. Muller et al, Phys. Rev. B (submitted). 
[19] T. Miyazaki et al, J. Phys. Soc. Jpn. 66, 2580 (1997). 
[20] N. Fujiwara et al., Phys. Rev. Lett. 80, 604 (1998). 



(b) 




J 

\ 



\ 




FIG. 1. (a) Schematic structure of the two-chain ladder 
compound SrCu203, which has Cu 2+ spin-1/2 magnetic mo- 
ments and antiferromagnetic exchange couplings J and J± 
along the chains and across the rungs, respectively. The 
isotropic inter-ladder coupling, indicated by the dashed lines, 
is frustrated; at site 2, for example, the coupling 1 - - 2 is can- 
celled by the coupling 2 - - 3. J c is the Cu-Sr-Cu coupling in 
the stacking direction of the ladders. For Sr(Cui_ z Zn z )203, 
a fraction z of the copper spins is randomly replaced by 
non-magnetic Zn 2 + ions, (b) Schematic of a single diluted 
ladder as investigated in the present Monte Carlo study. 
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FIG. 2. Correlation length versus degree of dilution of 
randomly diluted ladders. The temperatures were chosen such 
that T/J± = 3/100. Inset: Evidence for scaling of £Jj_/J in 
the variable zJ/J±. 
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FIG. 3. Correlation length and staggered susceptibility (in- 
set) for a randomly diluted ladder with Jj_/J = 1/2. 
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FIG. 4. Neel temperature T N (z) of Sr(Cui_ z Zn z ) 2 3 @ 
and points of constant correlation length £ = 18(1) of a single 
diluted ladder with J±/J = 1/2, as obtained from Monte 
Carlo. In this comparison, we use J = 1900-ftT JL6|. 
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